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(U)  ABSTRACT 


The  standard  Model  of  floating  probes  in  plasmas  is  reviewed. 
It  is  shown  that  a  Maxi  anna  charge  density  exists  for  the  applies* 
don  of  this  model.  An  explicit  expression  of  voltage  applied  to 
the  probes  in  terms  of  circuit  current  is  obtained.  Vhen  certain 
simplifying  assumptions  are  used,  this  equation  is  reduced  to  the 
Langmuir  single  probe  equation  and  the  Johnaon-Malter  double 
probe  equation. 

From  measurements  of  probe  current  versus  applied  probe  volt¬ 
age,  electron  temperature  at  each  respective  probe  can  be  deter 
mined.  One  example  is  given  to  illustrate  the  calculation. 
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IN  TRODU  CTION 


The  analysis  of  gas  discharge  processes  embodies  a  wide  variety 
of  theoretical  explanations  for  a  host  of  experimental  measurements. 
However,  present  day  knowledge  of  gas  discharge  properties  depends 
primarily  upon  experimental  measurements.  Perhaps  the  most  im¬ 
portant  variables  describing  a  gas  discharge  are  the  electron  and  ion 
densities,  and  their  velocity  distribution  functions.  Meaningful  meas¬ 
urements  of  temperature  of  plasma  constituents  require  them  to  be  in 
thermal  equilibrium.  When  this  condition  is  satisfied,  knowledge  of 
electron  and  ion  properties  as  determined  by  probes  is  required  to 
calculate  important  plasma  properties  such  as  ambipolar  diffusion 
coefficients,  excitation  rates,  and  ionization  rates.  There  are  at 
present  various  methods  for  analyzing  probe  data  which  are  valid  for 
special  experimental  sets  of  conditions.  Workers  in  this  area  of 
plasma  analysis  occasionally  invalidate  their  results  by  using  a  theory 
in  experimental  situations  where  the  assumptions  are  inapplicable.  It 
is  the  purpose  of  this  paper  to  show  the  connection  of  these  various 
probe  theories,  and  to  produce  a  generalized  theory.  A  simple  deriv¬ 
ation  of  probe  theory  is  presented  which  holds  for  various  probe  area 
ratios;  the  probes  are  not  required  to  be  surrounded  by  identical 
plasmas. 
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GENERAL  PROBE  THEORY 

Standard  treatment  of  experimental  data  taken  by  probe  methods 
allows  the  assumption  of  flat  plate  geometry  according  to  a  proof  by 
Druyvesten  (ref.  1),  who  concluded  that  any  convex  shape  for  a  probe 
should  be  appropriate.  The  plasma  sheath  condition  requires,  accord¬ 
ing  to  Langmuir  in  reference  2,  that  only  a  few  percent  of  electrons 
falling  on  the  sheath  may  collide  with  sheath  constituents  before  being 
collected  by  the  probe.  Hard  sphere  collisions  of  electrons  with  neutral 
particles  are  assumed  to  dominate  any  effect  of  scattering  interactions 
between  charged  particles  within  the  plasma.  The  limitations  on  this 
latter  assumption  will  be  discussed  in  appendix  A.  A  Maxwell -Boltzmann 
distribution  of  particle  speeds  is  assumed  so  that  the  concept  of  tem¬ 
perature  may  be  used.  Probes  are  considered  cold,  and  consequently 
do  not  emit  electrons. 

The  rate  of  sampling  charged  particles  by  the  probes  is  assumed 
to  be  small  compared  to  the  rate  of  charge  generation  of  the  local 
plasma,  so  that  the  plasma's  properties  outside  the  probe  sheath  will 
not  be  modified.  If  the  plasma  is  generated  from  thermal  excitation 
in  a  dense  gas,  the  conditions  required  for  this  assumption  are  possible 
to  obtain.  The  primary  condition  required  of  the  plasma,  (discussed 
in  more  detail  in  appendix  A),  is  that  the  plasma  be  diluted  by  neutral 
gas  atoms  in  thermal  equilibrium.  The  plasma,  consisting  only  of 
positive  ions  of  the  gas  species  and  free  electrons,  constitutes  a  very 
small  percentage  of  the  total  number  of  particles  present. 

The  basic  equations  governing  the  potentials  on  each  of  the  two 
probes  in  this  model  are  determined  from  Kirchhoff  s  law.  As  shown 
in  figure  1,  the  ion  currents  to  the  probes  are  ii  X  A  where  ii  is  the 
random  ion  current  in  the  plasma  around  the  probe  and  A  is  the  probe 
area.  A  similar  notation  for  random  electron  current  is  used,  in  which 
ig  is  the  random  electron  current.  Conventional  current  in  the  external 
circuit  is  denoted  by  I.  Primed  notations  are  used  to  refer  to  the  part 
of  the  total  random  currents  which  are  collected  by  a  probe  as  deter¬ 
mined  by  the  probe's  potential,  while  subscripts  refer  to  probe  1  and  2. 
Figure  1  shows  from  Kirchhoff' s  law  that 

iijAj  —  Ajiei  +  I  =  0,  (1) 

and 

iiaAj  “  Ajigj  ~  I  =  0,  (2) 
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where  Ai  and  A2  are  the  areas  of  thet/wo  probes.  Note  that  only  one 

I  >oo 

primed  term  per  equation  is  allowed  sittss«  ii  =•/,  f(v)dv  and 

.»  /  -v 

ie  ='co  f(v)dv  and  each  probe  must  bcttat  a  single  potential.  The 

function  f(v)  is  the  distribution  functioifl  for  x  component  of  velocity, 
where  V  =  mv*(%€<),  m  and  c  are  thesl'Uectron  mass  and  charge,  and 
v  is  the  velocity  of  the  electron. 

The  solutions  for  this  integration,  ili<nen  a  Maxwellian  distribution 
function  is  assumed  are 


Probe  A  is  assumed  to  be  mot  8  negative  than  probe  1 


Figure  1.  -  Schematic  diagram  n  of  probe  circuit. 
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exp  (-  ( V/kT^) 


(4) 


where  M  is  the  ion  mass,  T  the  temperature  in  degrees  Kelvin  and 
k  is  the  Boltzmann  constant.  Equation  (3)  may  be  substituted  into  (1) 
and  (2),  where  both  probes  are  assumed  to  be  at  negative  potentials 
with  respect  to  the  plasma  to  give  the  negative  potential  equations: 

Vi  =  Tei<k/«)  In  [(1  +  i/iuAiHiij/iei)].  (5) 

and 

V2  =  Te2<k/«)  In  1(1  -  1/ iij A2 ) ( ii2/ie2 )] .  (6) 

The  random  electron  current  terms  are  a  convenient  grouping  of 
constants  from  the  integration.  By  setting  I  equal  to  zero,  the  floating 
potential  Vf  for  any  surface  is  found  to  be 

Vf  =  l/2Te(kA)ln^g|iJ,  (7) 

which  is  in  its  standard  form,  assuming  thermal  equilibrium.  This 
potential  is  negative  with  respect  to  the  plasma  potential  when  |ie|  >  |ij| 
which  is  true  for  most  cases. 

Companion  equations  to  (5)  and  (6)  are  developed  when  the  probes 
are  assumed  to  be  at  some  positive  potential  with  respect  to  the  floating 
plasma  potential.  Equation  (4)  is  substituted  in  (1)  and  (2)  to  give  the 
positive  potential  equations 

V.=-T,<V.|ta  [(^)]  •  W 

and 

v,.-T«<k/.,i»[(i-ii)(ia.)]  .  (9> 

Here,  the  random  ion  current  terms  represent  a  convenient  grouping  of 
constants  from  the  integration.  The  positive  and  negative  potential 
equations  may  be  added  to  yield  complete  expressions  for  probe  poten¬ 
tials.  Measurements  of  circuit  voltage  represented  by  Va  in  figure  1, 
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correspond  in  an  experiment  to  differences  of  probe  potential,  as 
illustrated  in  figure  2.  Differences  between  the  complete  expressions 
for  probe  potentials,  taken  as  indicated  in  figure  2,  will  yield  the 
explicit  equation  of  Va  for  the  current-voltage  probe  characteristic: 


va  -  Tfcft/O  in  [(^)fe'))  -  T.W.I  ^ 

Tet(k/0  „  [(^)(^)]  -  Tei(k/«)  in  K^)(^)]  • 


(10) 


While  equation  (10)  holds  for  all  values  of  external  circuit  current  up 
to  its  limiting  value  as  determined  by  the  probe  area  and  Kirchhoff's  law. 
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only  two  terms  appear  in  the  calculation  at  any  one  time.  When  a  given 
probe  is  positive,  only  its  contributing  term  with  ion  temperature  as 
a  coefficient  appears  and  the  other  term  vanishes.  For  a  negative  probe, 
its  contributing  term  with  the  electron  temperature  as  a  coefficient 
appears  and  its  other  term  vanishes. 


MATHEMATICAL  CHARACTERISTICS  OF  THE  PROBE 
POTENTIAL  EQUATIONS 

It  is  important  to  recognize  the  special  properties  of  the  negative 
potential  equations  (eqs.  5  and  6).  The  following  definitions  will  be 
used  in  this  section  to  simplify  analysis: 


B‘  *  iiiA, 

(ID 

and 

B*  *  ifcAz 

(12) 

Potentials  of  each  probe  will  be  normalized  by  using  the  electron 
temperature  as  a  unit  of  voltage: 

V 

n‘Te 

(13) 

*  2 

"  Hi 

(14) 

Equations  (5)  and  (6)  may  be  rewritten  when 
with  the  above  definitions: 

Te  *  Tj  at  each  probe 

’l  l  =  In  | 

(1  +  B»)(m/M)1/2J 

(15) 

and 

’ll  *  In 

[(1  -  gBxMm/M)^] 

(16) 

where  g  =  Bj/Bi 


Figure  3  is  a  plot  of  and  n2  *or  hydrogen  where  the  ion  mass 
is  1 . 0  amu.  It  should  be  noted  that  is  limited  to  each  value  of  g 
determined  by  the  maximum  allowed  circuit  current.  For  values  of  g 
above  about  45,  the  small  probe  is  allowed  to  acquire  positive  potentials. 
For  g  =  100,  it  is  seen  from  the  figure  that  ^  is  limited  to  +  0.  87  for 
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Normalized  current  B 


Figure  3.  -  Normalized  current  versus  normalized  voltage 

on  each  probe. 


r)2  approaching  The  Langmuir  theory  (discussed  later)  is  valid 

for  hydrogen  only  when  t|2  is  nearly  equal  to  -  3.  76.  Large  area 
ratios  are  often  required  for  a  probe  to  develop  significant  positive 
potentials  consistent  with  equation  (23)  when  Tei  =  Tez. 

It  is  often  helpful  to  visualize  the  current  versus  voltage  on  a 
probe  for  a  given  species  of  particle  (ion  or  electron)  under  ideal  con¬ 
ditions.  Figure  4  is  a  plot  of  the  current  variation  with  probe  potential. 
Due  to  the  large  ratio  of  ion  and  electron  masses,  a  scaling  factor  was 
introduced  to  illustrate  each  current  on  the  same  plot. 

The  point  of  intersection  is  the  floating  potential  for  the  probe  in  a 
hypothetical  plasma  having  ie  =  4ij  (which  of  course,  is  not  realistic 
for  any  known  ion  at  the  same  temperature  as  the  electrons  in  its 


Figure  4.  -  Current  variation  with  probe  potential 

y2 

(M  =  16m)  Vf  =  k/c  In  *. 
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locality! )  Potential  zones  of  attraction  of  particles  show  all  random 
current  collected  by  the  probe  of  unit  area.  It  is  assumed  that  the 
sheath  dimensions  are  fixed.  Zones  of  repulsion  between  the  probe 
and  the  charges  show  the  collection  current  variation  according  to 
equation  (3)  and  (4). 


REDUCTION  TO  LANGMUIR  THEORY 


Consider,  as  an  example,  the  behavior  of  equation  (10)  when 
probe  2  acquires  a  positive  potential.  The  argument  of  the  logarithm 
of  the  first  and  third  term  in  equation  (10)  can  be  shown  to  be  equal  to 
1  by  using  Kirchhoff’s  law,  and  by  ascertaining  terms  that  are  fractions 
of  the  random  currents.  A  voltmeter  would  then  measure 


Va  =  k/« 


{-  T“ ln  [fee)©]  -  T“ ln  l(l^-)(i?)]] .  <17' 


as  the  external  circuit  voltage  for  external  circuit  current  I.  It  is 
apparent  that  equation  (17)  has  an  operating  bound.  The  maximum 
value  Imax  =  -iijA|  is  substituted  in  equation  (17).  The  required 
condition  is  that  the  second  term  must  always  be  positive  if  probe  2  is 
to  be  positive,  or 


0<T„ln[(l  )(!*»)].  »«> 

The  argument  of  the  logarithm  in  equation  (18)  must  therefore  be  less 
than  unity.  This  condition  establishes  a  requirement  on  both  probe 
areas  and  local  random  ion  currents  for  a  probe  to  acquire  a  positive 
potential,  or 


ii^>^-i.  d9) 

leiAJ  lei 


If  each  side  of  equation  (19)  is  equal,  then  probe  2  approaches  zero 
potential  when  probe  1  approaches 

An  important  case  arises  when 

Jil^l  »*»-!.  (20) 

ifcAj  i*2 

This  is  the  condition  for  a  "floating”  single  probe.  Equation  (18)  is 
transformed  into 


9 


(21) 


Va  *  -,/>T„<k/.)  ta(^;)  -  Tulk/O  m  [(l^)(lf)] 

for  the  current- voltage  relationship,  since  I  is  limited  by  ie2A2  which 
is  much  less  than  i^At. 

A  companion  to  equation  (21)  is  developed,  when  the  NfloatingN  probe 
goes  negative  with  respect  to  its  surrounding  plasma,  by  a  similar 
treatment  of  equation  (10)  for  -  I: 

v.  *  -  V,Tel(k/. )  ln(2jlu.)  -  T«(k/«  ,  1«[(^)(^)]  (22) 

The  first  term  in  equations  (21)  and  (22)  is  recognized  to  be  the  floating 
probe  potential  for  probe  1.  By  making  use  of  this  fact  and  substituting 
in  the  Kirchhoff  relation  for  I,  equations  (21)  and  (22)  may  be  written 
in  a  different  form: 


V2  =  Va  +  Vfl  =  TijtkA)  In  {iu)  -  T^k/c)  In  (i!,) 
V2  =  Va  +  Vft  =  TeaOc/c)  In  (i*)  -  T^k/c)  In  (i^) 


(23) 


Equations  (23)  were  first  derived  by  Langmuir  and  Mott-Smith  (ref.  3) 
for  "floating"  single  probes.  Thus,  under  the  approximations  shown 
above,  equation  (10)  may  be  reduced  to  the  special  case  of  the  Langmuir 
single  probe  equation.  The  potential  of  the  floating  probe  with  respect 
to  the  surrounding  plasma  "zero"  potential  is  plotted  as  a  function  of 
the  natural  logarithm  of  the  current  falling  on  the  probe.  Since  the 
first  term  in  each  of  equations  (23)  is  a  constant,  the  straight  line 
slope  of  the  current-voltage  data  is  both  a  determination  of 
the  Maxwellian  character  of  the  collected  charges  and  their 
temperature. 


REDUCTION  TO  JOHNSON  AND  MALTER  THEORY 


With  equation  10  as  the  starting  point,  it  is  assumed  that  both 
probes  have  area  ratios  which  fail  to  satisfy  the  conditions  prescribed 
by  equations  (19)and  (20).  The  measured  voltage  between  the  probes 
would  then  depend  primarily  upon  the  temperature  and  random  current 
of  the  local  electrons  surrounding  each  respective  probe,  as  seen  in 
equation  (24). 


=  -Tel(k/«)  In 


+  Tej(k/«)  In 


(24) 
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For  convenience,  the  Johnson  and  Malter  model  (ref.  3)  assumed  the 
probes  were  so  close  together  that  no  significant  difference  between 
local  electron  temperatures  would  be  permitted.  Equation  (17)  is 
written  in  its  simplest  form  by  using  these  conditions. 

v,  =  T.(k/.)  <25> 

Differentiation  of  equation  (25)  results  in 


Te(k/«) 


_ (Aitji  +  Agiig) _ 

(A2ii2  -  l|  )(A,ii7  +  II 

Va  =  0  1  Va 


(26) 


where  R0  was  defined  as  the  "equivalent  resistance"  of  the  plasma. 
The  following  substitutions  are  required  to  place  equation  (26)  in  a 
familiar  form:  Lip  =  Aji^  +  Aji^  and  G  *  [A«l5s2/  *  iplv.  =0 

Electron  temperature  may  then  be  calculated  by 


Te  =  1 1 , 600(G  -  Gz)Rq£  £p 


(27) 


This  result  was  first  derived  by  Johnson  and  Malter  (ref.  3)  as  the 
"equivalent  resistance  method".  Use  of  equation  (27)  has  an  important 
restriction  which  is  worth  noting.  The  calculated  electron  temperature 
is  some  sort  of  effective  or  average  value  when  the  probes  are  placed 
in  a  plasma  of  unequal  local  electron  temperatures.  For  this  situation, 
Te  does  not  represent  the  actual  value  at  either  probe. 


GENERAL  DOUBLE  PROBE  METHOD 

The  expression  for  the  current* voltage  relation  of  equation  (10) 
may  be  also  used  for  the  general  case  in  which  Tej  V  T e2  and 
HiAi  ^  ii*Az.  The  equation  for  calculating  these  temperatures  from 
current- voltage  data  is  found  by  expanding  equation  (24)  to  reflect  the 
floating  potential  of  each  probe  according  to  equation  (7): 

Va  =  Tel(k/«)  In  -  Tel(k/«)  In  -  Vfl  +  Vfc.  (28) 

The  terms  Vf2  “  Vfi  are  the  differences  in  floating  potentials  between 
probes  1  and  2,  and  will  be  denoted  as  AVf  as  shown  in  figure  5. 
When  the  external  circuit  current  is  zero,  the  voltage  measured  is 
the  difference  in  floating  potentials.  In  order  to  calculate  T^  and 
Tez  from  actual  data,  a  knowledge  of  iiA  for  each  probe  is  required. 
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Figure  5.  -  Typical  double  probe  data  characteristic  curve. 

This  quantity,  ijA,  is  the  saturated  current  and  should  be  almost 
constant  for  all  values  of  Va. 

The  following  table  is  helpful  in  identifying  which  saturated  cur¬ 
rent  refers  to  which  probe,  and  the  appropriate  signs  of  current- 
voltage  for  analysis  of  data. 

An  analysis  of  data  with  unequal  electron  temperatures  is  given 
in  appendix  Bto  illustrate  the  procedure. 


TABLE  I 


TABLE  OF  SATURATED  CURRENT  IN  REFERENCE 
TO  PROBE  FOR  DATA  ANALYSIS 


Conditions 

Tei  >  Te2 

Tet  >  Tel 

H 

II 

o 

Va  is  positive 

Va  is  negative 

va  -  +  00 

I  -  -  iuA, 

I  -  -  i^Ai 

va  -  -  - 

I  -•>  i^Aj 

I  -*  +  ifcAj 

va  *  0 

H 

II 

+ 

H 

II 

o 

1  =  "  1  ^va=o| 

CONCLUSION 


A  general  probe  theory  has  been  derived  which  may  be  applied  to 
data  collected  from  floating  probes  (see  eq.  10)  in  dilute  plasmas. 
This  theory  includes  the  effects  of  the  probes  being  placed  in  plasma 
regions  of  unlike  temperature  environments.  From  an  explicit 
expression  of  measured  voltage  in  terms  of  measured  circuit  current, 
it  is  shown  that  with  appropriate  simplifications  the  Langmuir  "single 
probe"  equation  and  the  Johnson-Malter  equation  are  obtained. 
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Appendix  A 

TEMPERATURE  LIMITATION  OF  DOUBLE  PROBE  THEORY 

A  theory  for  determining  the  temperature  of  a  plasma  by  electri¬ 
cal  probes  is  described  in  the  main  section  of  this  report.  It  is  valid 
for  electrons  and  ions  which  are  scattered  by  neutral  atoms  only,  since 
no  consideration  was  made  for  Coulomb  interaction  between  the 
charges.  It  becomes  important  to  determine  the  temperature  for 
which  Coulomb  forces  between  charges  become  significant  in  modifying 
the  theory.  Consideration  of  the  electrical  conductivity  as  a  function 
of  temperature  will  illustrate  the  maximum  allowable  temperature  for 
which  confidence  may  be  placed  in  the  interpretation  of  probe 
measurements  by  the  above  mentioned  theory. 

Calculation  of  electrical  conductivity  in  a  gas  requires  knowledge 
of  the  mobility  of  the  change  carriers.  The  mobilities  of  charged 
particles  in  a  partially  ionised  gas  are  influenced  both  by  Coulomb  and 
collision  forces.  Electrical  conductivity  maybe  easily  calculated  for 
two  extreme  cases,  that  of  very  little  ionization  and  the  other  of 
almost  complete  ionisation.  The  isobaric  conductivities  calculated 
for  both  cases,  assuming  thermal  ionization,  increase  monotonically 
with  temperature.  Experimental  results  agree  very  well  with  the 
theory  for  each  case.  For  the  intermediate  degree  of  ionization,  there 
'is  a  transition  zone  in  which  mobilities  are  influenced  by  both  types  of' 
collision  forces.  The  effluent  gases  from  a  plasma  jet  are  of  the 
proper  temperature  and  density  ranges  to  manifest  conductivities 
characteristic  of  all  three  regions:  collision,  transition,  and 
Coulomb  forces. 

An  expression  for  the  electrical  conductivity  of  a  gas  in  the  low 
temperature  region  (valid  to  7,000*K)  was  given  by  Lin,  (ref.  4)  as 
follows: 

<r=nee*r/me,  (29) 

where  v  is  the  gas  conductivity,  t  the  mean  collision  time  for 
electron-atom  (or  molecule)  collisions,  and  ne,  e,  and  me  are  the 
electron  density  (no.  /cm1 2)  the  electron  charge,  and  the  electron  mass, 
respectively.  The  following  assumptions  are  implicit  in  the  derivation 
of  equation  29: 

1.  In  the  region  considered,  the  gas  density  is  large  enough  to 
provide  sufficiently  numerous  collisions  between  charged  particles 
and  neutrals  so  that  statistical  equilibrium  is  attained; 

2.  No  electron  attachment  or  other  long-lived  interactions  occur 
between  neutral  and  charged  particles; 
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3.  Values  of  the  temperature  and  electron  density  are  such  that 
the  Debye  shielding  distance  is  small  in  comparison  to  any  critical 
physical  length  such  as  probe  spacing  or  size; 

4.  Electrical  neutrality  exists; 

5.  No  external  electric  or  magnetic  fields  exists;  and 

6.  Due  to  high  electron  mobility  relative  to  the  ions,  the 
electrons  are  the  current  carriers. 


By  appropriate  evaluation  of  t  from  kinetic  theory,  equation  (29) 
may  be  written  for  argon: 

<r  =  5.  4  X  107 — ^-p7(mho/meter)  (30) 

ntTeVz 


where  nt  is  the  number  density  of  neutral  particles,  and  Te  is  the 
electron  temperature.  The  atomic  radius  of  argon  (given  in  the 
Handbook  of  Chem.  and  Phys.)  was  used  in  evaluating  the 
proportionality  factor,  5.4  X  10T,  in  equation  (30)  and  ne  is 
determined  by  Saha's  ionization  equation.  For  the  case  cited 


ne  =  4.66  X  10< 


|b'(T) 

b(T) 


1. 82  X  10* 
T 


.electron. 
(  rJ  ) 


(31) 


where  b'(T)  and  b(T)  are  the  partition  functions  for  singly  ionized 
and  neutral  argon,  respectively,  (ref.  5).  By  assuming  the  perfect 
gas  state,  the  isobaric  conductivity  for  "dilute"  plasma  may  be 
determined  as  a  function  of  temperature.  Figure  A-l  shows  graph¬ 
ically  the  results  of  such  a  calculation  for  the  1  cm  Hg  isobar  in  an 
argon  plasma. 


For  fully  ionized  or  "dense"  plasmas  in  which  Coulomb  collisions 
are  expected  to  predominate,  an  expression  for  plasma  resistivity 
has  been  developed  by  Spitzer  and  HSrm  (ref.  6) 


R  = 


3.8  X  10*ZlnA 
Ts/z 


(ohm-meter) 


(32) 


where  Z  is  the  ion  charge  in  multiples  of  the  electron  charge  and 
In  A  is  a  function  which  is  tabulated  by  Spitzer  and  HSrm  (ref.  6). 

The  values  of  plasma  conductivity  calculated  from  equation  (32)  are 
shown  in  figure  A-l,  together  with  those  values  calculated  from 
equation  (30).  It  may  be  seen  from  this  graph  that  for  an  argon  plasma 
at  a  pressure  of  1  cm  Hg,  the  probe  theory  may  be  applied  for 
temperatures  up  to  the  transition  zone  from  "dilute"  to  dense  which 
is  6,  000*K  to  8,  500*K. 
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Appendix-  B 

AN  ANALYSIS  OF  SOME  PROBE  DATA  BY  THE  GENERAL 
DOUBLE  PROBE  THEORY 


Johnson  and  Malter  (ref.  3)  have  obtained  double  probe  data  in 
quiescent  decaying  plasmas.  Their  data  were  obtained  from  plasmas 
generated  by  a  pulse  discharge.  Probe  currents  versus  time  were 
measured  at  fixe<?  values  of  potential  difference  between  probes.  The 
data  shown  in  figure  A-l  are  taken  directly  from  their  paper. The  points 
shown  on  the  curve  1  were  taken  400  sec  after  the  plasma  generating 
discharge  had  closed.  This  period  was  determined  by  Johnson  and 
Malter  to  be  sufficient  to  establish  equilibrium.  The  area  of  each 
probe  in  their  experiments  was  equal  so  that  the  ratio  At/A2  which 
appears  in  the  simplified  probe  theory  is  unity  for  this  discussion. 

For  Aj/Aj  less  than  a  critical  value  which  is  determined  by  the 
ion  mass  in  the  plasma  (see  equation  20),  the  general  probe  theory 
gives  the  applied  voltage -cur rent  relation  by  the  equation 

Va  =  k/«Teiln(l-l/iilA1)-k/€Te2ln(l  +  l/ii2A2)  +  AV£  (33) 

where  AVf  is  the  difference  in  plasma-probe  potential  between  the 
two  probes. 

In  order  to  use  the  data  of  Johnson  and  Malter  in  the  above 
equation  constant  values  of  i|A  must  be  assumed.  Observations  of 
most  data  reveals  the  continued  increase  in  I  after  the  supposed 
saturation  value  of  I  has  been  reached.  This  increase  in  I  is 
due  to  expansion  of  the  sheath  with  voltage.  The  continued  increase 
in  I  can  be  attributed  to  an  effective  increase  of  probe  area,  which 
occurs  because  of  an  increase  in  sheath  volume  produced  by  an  in¬ 
crease  in  current  to  the  probe.  For  a  first  approximation,  this  per¬ 
turbation  is  assumed  linear  with  I.  The  dashed  lines  on  figure  B-l 
show  how  the  magnitude  of  the  perturbation  has  been  estimated.  The. 
figure  also  shows  the  data  of  Johnson  and  Malter  as  modified  by  the 
perturbation  correction.  The  modified  data  points  were  obtained  by 
reducing  the  current  I  by  an  amount  equal  to  the  perturbation.  For 
example,  in  figure  B-l,  the  estimated  perturbation  on  I  at  data  point 
£  is  d*  as  shown  by  the  solid  bar  between  the  dashed  lines.  There.- 
fore,  point  a  is  lowered  by  an  amount,  d,  and  placed  on  the  mod¬ 
ified  data  chart  of  figure  B-2.  A  similar  correction  was  made  for  17 
other  points. 

The  modified  Johnson  and  Malter  data  were  analyzed  by  the  least 
squares  method  to  fit  the  equation 
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yi  =  ai  +  a2ln(l  -  Xi/xj)  +  a3ln(l  +  xi/x2) 


(34) 


which  is  of  the  same  form  as  equation  (32).  The  constants  a2,  a2, 

and  a3  are  the  unknowns  to  be  determined.  The  x^  and  are  the 
given  data  points.  The  values  of  X]  and  x2  are  determined  graphi¬ 
cally  (see  fig.  B-2). 

The  equations  which  explicitly  satisfy  the  least  square  minimiza 
tion  requirement  for  the  problem  cited  here  ares 


^  =  m  [a*  + 

0=  S  |a2  +  a2ln(l 

0=2  [a,  +  a2ln(l 

m  l 


-  xi/xj)  +  a3ln(l  +  xi/x2) 

-  x-/ x^)  +  a2ln(l  +  xi/x2) 

-  Xi/x^  +  a3ln(  1  +  Xi/x2) 


Yi]  (35) 

yj  ln(l  -  x j/ Xj )  (36) 
y[]  ln(l  +  Xi/x2)  (37) 


Eighteen  points  were  taken  from  the  modified  Johnson  and  Malter 
data.  These  points  were  chosen  more  or  less  uniformly  along  the 
current -voltage  curve.  The  value  of  m  in  equations  (35)  through  (37) 
therefore,  runs  from  1  through  18.  It  may  be  noted  here  that  the 
present  analysis  does  not  fully  utilize  the  data  available,  since  there 
were  30  available  data  points.  The  values  of  the  coefficients  obtained 
from  the  simultaneous  solution  of  equations  (35)  through  (37)  are: 
a,  8  •  0.  0752d  volt,  aa  =  0.  07253  volt,  and  a3  =  -  0.10045  volt.  The 
graphically  determined  values  of  xj  and  x2  were  used  in  equations 
(35)  through  (37);  these  values  are:  x2  =  0.3723  amp  and  x2  =  0.6012 
amp.  By  comparing  the  coefficients  aj,  a2,  and  a3  with  their  corre¬ 
sponding  terms  in  equation  (32),  the  electron  temperatures  at  the 
probes  and  the  difference  in  plasma-probe  potentials  were  determined 
to  be  Te  =  841'K,  Tea  =  1,  165°K,  and  AVf  *  -  0. 07528  volt.  The 
curve  in  figure  B-2  is  generated  from  equation  (34)  after  appropriate 
conversion  of  coefficients  a2,  a2,  and  a3. 


From  the  Langmuir  theory  as  modified  and  used  by  Johnson  and 
Malter,  (ref.  3)  the  values  of  Te{  and  TQ  are  not  obtained  separately. 
However,  they  obtained  an  average  value  o/  Te,  the  electron  temper¬ 
ature,  of  Te  =  1,015°K.  This  value  compares  favorably  with  the 

arithmetic  average  of  T_  and  T  which  is  1 , 003*K. 

ei  e2 
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Dr.  J.  L.  Martin 

Dr.  Aba  M.  Zarsm 
Electro -Optical  Syatsma 
1 70  N.  Daiay  Avanua 
Pasadena,  California 

Dr.  William  Clohssay,  Diractor 

Technical  Davalopmsnt 

Tha  Martin  Company,  Danvar  Division 

P.  O.  Box  179 

Danvar  1,  Colorado 

Mr.  rrank  Bradley,  Diractor 
Requirements 

Tha  Martin  Company,  Baltimors  Division 
Baltimore  3,  Maryland 

Radio  Corporation  of  America 
Mooreetown,  New  Jersey 
ATTN:  Mr.  Alvin  C.  Gottlieb 

Rocketdyne 

6633  Canoga  Avanua 

Canoga  Park,  California 

ATTN;  Dr.  Robert  J.  Thompson,  Jr. 

Dr.  C.  J.  McDole 


Vidya 

2626  Hanover  Street 
Palo  Alto,  California 
ATTN:  Dr.  M.  G  Boobar 

Saint  Louie  Univaraity 

221  N.  Grand,  St.  Louie,  Missouri 

ATTN:  Mr.  A.  H.  Weber,  Phyaica  Dept. 

Project  MATTERHORN 
Jameo  For  rental  Research  Canter 
Princeton  University,  P.  O.  Box  451 
Princeton,  New  Jersey 

Dr.  C.  W.  Tombaugh 
Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture  and  Mechanical  Arts 
Box  5 41,  State  College,  New  Mexico 

Professor  Kart  Wohl 
Chemical  Engineering  Department 
University  of  Delaware 
Newark,  Delaware 

University  of  Alabama 

University,  Alabama 

ATTN:  Dr.  F.  H.  Mitchell,  Boa  1392 

Bureau  of  Engineering,  bm  6162 
Professor  H.  Kuan  sol.  Department  of 
Mechanical  Engineering 
Department  of  Chemical  Engineering 

Dr.  Loren  E.  Bollinger 
Rocket  Research  Laboratory 
Department  of  Aeronautical  Engineering 
The  Ohio  State  University 
Columbus  10,  Ohio 

Dr.  Nter  Debye 
Chemistry  Department 
Cornell  University 
Ithaca,  New  York 
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Dr.  Paul  Groaa 
Department  of  Chemietry 
Duke  University 
Durham,  North  Carolina 

Dr.  J.  Richard  Haskins 
Gettysburg  College 
Gettysburg,  Pennsylvania 

Dr.  William  D.  Will  lama 
Harding  Collage 
Box  602 

Searcy,  Arkansas 

Prof.  Paul  W.  Gills s 
Department  of  Chemietry 
The  University  of  Kansas 
Lawrence,  Kansas 

Dr.  Harold  G.  Donnelly,  Head 

Chemical  and  Metallurgical  ftigineering  Department 
Wayne  State  University 
Detroit  2.  Michigan 

Armour  Research  Foundation 
Ulinoie  Institute  of  Technology 
Technical  Canter 
10  W.  35th  Street 
Chicage  16,  Ulinoie 

Rohm  and  Haas  Company 

Rada  tone  Arsenal  Research  Division 

Redstone  Arsenal,  Alabama 

ATTN:  Librarian 

Thiokol  Chemical  Corporation 
Redstone  Division 
Redstone  Arsenal,  Alabama 
ATTN:  Librarian 

George  C.  Marshall  Space  Flight  Center 
National  Aeronautics  and  Space  Administration 
Redstone  Arsenal,  Alabama 
ATTN:  M-SAT-DIR  (Dr.  O.  H.  Lange) 

M-S  and  M-M  (Dr.  W.  R.  Lucas) 

M-RP-T  (Mr.  G.  B.  HeUer, 

Mr.  D.  W.  Gates) 

National  Aeronautics  and  Space  Administration 
Manned  Spacecraft  Center 
Langley  AFB,  Virginia 
ATTN:  Mr  R.  B.  Erb 

Mr.  J.  Kotaachik 


Commanding  General 
U.  S.  Army  Ordnance  Missile  Command 
Redstone  Arsenal,  Alabama 
ATTN:  ORDXM-XS 

OR0XM-XE 
ORDXM-XGR 
ORDXM-R 
ORDXM-RE 

OADXM-REX  (Mr.  E.  T.  Carr) 

ORDXM-RF 

ORDXM-RC 

ORDtXM-RP 

ORDXM-RH 

OftDXM-Rl 

OKDXM-RM 

OADXM-RK 

ORDXM-RS 

ORDXM-RG 

ORDXM-RR 

ORDXM-RB 

OADXM-RAP 
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